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Aggregation of two oppositely charged colloids in salty water is studied. We focus on the role 
of Coulomb interaction in strongly asymmetric systems in which the charge and size of one colloid 
is much larger than the other one. In the solution, each large colloid (macroion) attracts certain 
number of oppositely charged small colloids (Z-ion) to form a complex. If the concentration ratio 
of the two colloids is such that complexes are not strongly charged, they condense in a macroscopic 
aggregate. As a result, the phase diagram in a plane of concentrations of two colloids consists of 
an aggregation domain sandwiched between two domains of stable solutions of complexes. The 
aggregation domain has a central part of total aggregation and two wings corresponding to partial 
aggregation. A quantitative theory of the phase diagram in the presence of monovalent salt is 
developed. It is shown that as the Debye-Hiickel screening radius r s decreases, the aggregation 
domain grows, but the relative size of the partial aggregation domains becomes much smaller. As 
an important application of the theory, we consider solutions of long double-helix DNA with strongly 
charged positive spheres (artificial chromatin). We also consider implications of our theory for in 
vitro experiments with the natural chromatin. Finally, the effect of different shapes of macroions 
on the phase diagram is discussed. 

PACS numbers: 61.25.Hq, 82.70.Dd, 87.14.Gg, 87.15. Nn 



I. INTRODUCTION 

Aggregation or self-assembly of oppositely charged col- 
loids is a general phenomenon in biology, pharmacol- 
ogy and chemical engineering. The most famous bi- 
ological example is the chromatin made of long nega- 
tively charged double-helix DNA and positively charged 
histone octamcrs 1 . Due to Coulomb interaction, DNA 
winds around many octamcrs to form a beads-on-a-string 
structure, also called 10 nm fiber (see Fig.^for an illus- 
tration). This 10 nm fiber may self-assemble into a 30 nm 
fiber, which is the major building material of a chromo- 
some. The formation of 30 nm fiber strongly depends on 
the concentration of salt in the solution. This means that 
the Coulomb interaction plays a crucial role 2 . The best 
known pharmacological example is the problem of gene 
therapy. In this case, a negatively charged DNA helix 
should penetrate through a negatively charged cell mem- 
brane. To do this, DNA has to be neutralized or over- 
charged by complexation with positive polyelectrolytes 
or colloids. At the same time aggregation of these com- 
plexes can be useful or should be avoided2^iSiL2iS. In- 
dustrial examples of aggregation include using cationic 
polyclcctrolytc as coagulants for paper manufacturing, 
mineral separation, and the aggregation-induced removal 
of particulate matter from the aqueous phase in water 
and waste water treatment processes^. Given such wide 
applications, it is therefore interesting to construct a gen- 
eral physical theory on aggregation of oppositely charged 
colloids. 

In this paper, we consider the equilibrium state of 
two kinds of oppositely charged colloids in salty wa- 
ter and construct the phase diagram of such a sys- 
tem. Without losing generality, we call the larger col- 
loid "macroion" with negative charge —Q, concentration 




FIG. 1: A beads-on-a-string structure of the complex of a 
long negative polymer macroion with positive spherical Z- 
ions (artificial chromatin). 



p, and the smaller colloid "Z-ion" with positive charge 
Ze (e is the proton charge), concentration s. Macroions 
can be big spheres (big colloid particles), rigid cylinders 
(short DNA) or long semi-flexible polymers (long DNA). 
Z-ions can be small spheres (nucleosome core particles, 
micelles or dendrimers) or short polymers (polyamines). 
Most of such systems are strongly asymmetric in the 
sense that the size and charge of the macroion are much 
larger than those of the Z-ion (Figs. ^ 12 GJ>- Therefore, 
defining the number of Z-ions neutralizing one macroion 
as Ni = Q/Ze (the subscript i denoting "isoelectric"), 
we focus on systems with Ni 3> 1. 

In the previous paper— this problem was considered 
for long semi-flexible polymers (DNA double helices) as 
macroions and rigid synthetic spheres with very large 
charge as Z-ions (Fig. ^| without monovalent salt. In 
such a system each polymer macroion winds around a 
number of spherical Z-ions and forms a periodic necklace- 
like structure which is similar to the natural chromatin. 
We call it "artificial chromatin". The phase diagram 
Fig. 0]} was obtained in a plane of the macroion (DNA) 
concentration, p, and the Z-ion (spheres) concentration, 
s (notations s and p are introduced as abbreviation for 
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FIG. 2: A complex of a negative spherical macroion with 
positive spherical Z-ions condensing on it. 




FIG. 3: A complex of a negative cylindrical macroion with 
positive spherical Z-ions condensing on it. Z-ions form 2D 
Wigner-crystal-like liquid on the surface of the macroion. 



"sphere" and "polymer"). As seen from Fig.Bi around 
the line of neutrality (the dashed line) where macroion- 
Z-ions complexes are almost neutral, they condense in a 
macroscopic aggregate (the gray region). Far away from 
this line, the complexes are strongly charged and they 
stay free in the solution (the white region). 

It should be emphasized that strong correlations play 
a crucial role in the origin of this picture. First, spherical 
Z-ions in the necklace repel each other to form a ID liquid 
with almost periodic structure similar to one-dimensional 
Wigner crystal (see Fig. . This leads to a correlation 
energy gain in addition to the mean field Coulomb en- 
ergy of the spheres. The difference between this gain 
and the entropic gain in the bulk per Z-ion plays the 
role of a voltage which may overcharge necklaces (mak- 
ing them positive). As a result, a free macroion- Z-ions 
complex can be either positive (above the dashed line in 
Fig- BO or negative (below the dashed line in Fig. BO de- 
pending on concentrations s and p. Second, the polymer 
turns wound around a spherical Z-ion repel each other 
and form an almost equidistant coil. This correlation 
offers a short-range attraction between complexes as il- 
lustrated by Fig. [5] When complexes are almost neutral, 
this attraction is able to condense them. When they are 
strongly charged, Coulomb repulsion is larger than the 
short-range attraction and they stay free. 

An important feature of Fig. B is two relatively wide 
domains, where aggregation of macroion- Z-ions com- 
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FIG. 4: Phase diagrams of the artificial chromatin system, a: 
without monovalent salt, b: with monovalent salt. (Spherical 
and cylindrical macroion systems have similar phase diagrams 
with a much larger Sd-) p is the concentration of macroions 
(DNA), s is the concentration of Z-ions (spheres). Plus and 
minus are the signs of the charge of free macroion- Z-ions com- 
plexes. The dark gray region is the domain of total aggrega- 
tion of complexes, while the light gray region is the domain 
of their partial aggregation. The white region is the domain 
of free complexes. s c and Sd are concentrations of Z-ions at 
the boundary of the aggregation domain when p — ► 0. At 
the dotted "isoelectric" line the absolute values of the total 
charges of macroions and Z-ions in the solution are equal. At 
the dashed "neutrality" line free complexes are neutral. 



plexes is only partial (the light gray region) . This means 
that certain fraction of complexes aggregates, while oth- 
ers stay free in the solution. The reason for partial aggre- 
gation is redistribution of Z-ions between aggregated and 
free complexes. When (p, s) is such that the Coulomb 
repulsion between complexes is slightly larger than the 
short-range attraction, Z-ions can redistribute them- 
selves so that for a fraction of complexes the Coulomb 
repulsion becomes substantially smaller and they aggre- 
gate, while for the rest of complexes the Coulomb repul- 
sion becomes even larger and they stay fresi^. With- 
out monovalent salt, the macroscopic aggregate should 
be practically neutral (a net charge defined as the sum 
of charges of the macroion and Z-ions should be close 
to zero). As a result, the domain of total aggregation is 
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FIG. 5: Cross section through the centers of two touching 
positive spherical Z-ions with turns of negative semi-flexible 
polymers (gray) wound around them. At the place where 
two spheres touch each other (the rectangular), the density 
of polymer segments doubles and the correlation energy is 
gained. 



narrow while the domain of partial aggregation is much 
wider (see Fig. QJi). 

In this paper, we develop a general phenomenologi- 
cal theory of complexation and aggregation of oppositely 
charged macroions and Z-ions in the presence of mono- 
valent salt. In order to formulate this theory we use the 
simplest model of spherical rigid macroions and small 
Z-ions shown in Fig. [3 In the absence of monovalent 
salt, the phase diagram of this system looks like Fig. 0Ji, 
where p is the concentration of spherical macroions and 
s is the concentration of Z-ions. Screening by monova- 
lent salt effectively truncates the Coulomb interaction at 
the Debye-Huckel screening radius r s . If r s is smaller 
than the size of the spherical macroions, the Coulomb 
repulsion of two macroion- Z-ions complexes is almost 
completely screened even when they touch each other. 
Therefore, aggregation is possible even if macroion-Z- 
ions complexes carry a net charge (defined as the sum 
of charges of the macroion and Z-ions). As a result, 
screening strongly modifies the phase diagram as shown 
in Fig. 0Jd. Comparing with the phase diagram without 
monovalent salt (Fig. 0^), we see that the aggregation 
domain grows. At the same time, the relative size of the 
partial aggregation domain (the light gray region) in the 
whole aggregation domain decreases, because redistribu- 
tion of Z-ions between aggregated and free complexes 
becomes less important. 

Our phenomenological theory of Sec. II contains sev- 
eral parameters depending on charges, sizes, shapes and 
flexibility of macroions and Z-ions. In sections III, IV 
and V we evaluate microscopically these parameters for 
specific pairs of macroions and Z-ions. In Sec. Ill we 
do this for a system of spherical macroions and small 
Z-ions (Fig. |2J) and find slopes of phase diagram bound- 
aries in this case. In Sec. IV we repeat this calculation 
for cylindrical macroions (Fig. QJ. In Sec. IV we return 



to artificial chromatin (Fig.^l and revise some results of 
Ref. ^ related to the screening effect of monovalent salt. 

We find that although all these systems have phase dia- 
gram similar to Fig.^Ja, the width of aggregation domain 
and relative width of the partial aggregation domains de- 
pend on the pair. Namely, rigid cylindrical macroions 
have the largest aggregation domain, because parallel 
cylinders have relatively larger area where they can touch 
each other. Artificial chromatin (Fig. ^) on the other 
side has the smallest aggregation domain. This peculiar- 
ity will be explained in Sec. V. There we also show that 
phase diagram of Fig. seems to agree with results of 
experiments in vitro designed to understand regulation 
of natural chromatin^. 



II. PHENOMENOLOGICAL THEORY OF 
AGGREGATION IN THE PRESENCE OF 
MONOVALENT SALT 

In this section we discuss the general theory of the 
phase diagram in the presence of monovalent salt using 
the system of spherical macroions shown in Fig. [21 Ag- 
gregation of two such spherical macroions is illustrated 
by Figs. EJand Z-ions form a strong correlated liquid 
on the surface of the macroion. The negative correlation 
energy in this liquid results in the voltage which may 
overcharge the complex. It also induces a short-range 
attraction between complexes since in the spot where 
two complexes touch each other (see Fig. |SJ| the surface 
density of the correlated liquid is doubled and the en- 
ergy per Z-'vm is reduced 16 . This attractive force de- 
cays as exp(— 2-7T(i/v / 3V) (for a triangular lattice)^ with 
the distance d between two surfaces, where r is the half 
of average distance between nearest neighbor Z-ions on 
the macroion surface, or, in other words, the radius 
of the Wigner-Seitz cell. Consequently, only a small 
disk-like contact region with radius W — y/jRr, where 
7 = v3/27T ~ 0.3, contributes to the attraction of two 
complexes (see Fig. EJ). We call it "sticky region". 

Similarly to artificial chromatin, this correlation at- 
traction leads to phase diagram Fig. 2t m the absence 
of screening. We should remember, however, that in 
the case of Fig. s stands for the concentration of 
small spherical Z-ions, and p is the concentration of large 
spherical macroions. We show below that screening not 
only enlarges the aggregation domains, but also leads to 
relatively smaller size of the partial aggregation domain. 
These effects become strong when r s is smaller than R. 
They grow when r s is reduced all the way to r and even 
at r s < r. Only when r s is so small that interaction of a 
Z-ion with macroion surface becomes smaller than fcsT, 
most of Z-ions leave the macroion surface, and aggrega- 
tion becomes impossible. Below we concentrate on the 
case r <C r s <C R which can be described by simple and 
universal phenomenological theory. 

There is a key difference between the cases without 
screening discussed in Ref. and with screening effect 
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FIG. 6: The cross section of the touching area of two aggre- 
gated spherical macroions. The sticky region is shown by the 
rectangular. In this region, the density of Z-ions doubles and 
the correlation energy is gained. 



of monovalent salt. In the first case, each macroion of 
the aggregate can only carry Ni Z-ions, i.e., the aggre- 
gate is almost neutral. This is because if macroion- Z- 
ions complexes were charged, the Coulomb potential of 
the aggregate would increase with its size and become 
too large to continue aggregation. On the other hand, 
in the presence of strong screening by monovalent salt, 
when the local distance between the two touching sur- 
faces, d, is much larger than r s , the Coulomb interaction 
between two surfaces is exponentially small and can be 
neglected (see Fig. EJ . It only affects the disk-like surface 
region with radius W s = \/2Rr s located inside the cylin- 
der where d < 2r s . We refer to such a cylinder as "r s - 
cylinder" (shown by the dashed rectangular in Fig. [7J. 
Therefore, if we count only charges of macroions and Z- 
ions, the aggregate of macroion- Z-ions complexes need 
not be neutral. Most of the surface area of these com- 
plexes can be overcharged or undercharged in the same 
way as surfaces of free complexes. 

We argue that macroion surface charge density, a(d) 
(including charge of the macroion and Z-ions), inside r a - 
cylinder is much smaller than the one outside r s -cylinder, 
(Tq- In electrostatics, it is known that for two charged con- 
tacting metallic spheres in vacuum, the surface charge 
density goes to zero when we move to the point of con- 
tact (one can show that it decays as exp(—-Ky/R/d)). Our 
case is similar to this example in the sense that Z-ions 
equilibrate on macroion surfaces so that everywhere they 
acquire the same electrochemical potential, p c {N) + e<f>. 
The chemical potential fJ. c (N) determined by correlations 
of Z-ions depends on concentration of Z-ions on the sur- 
face of macroion, which, as shown below, can be regarded 
as a constant. This means that the electrostatic potential 
4> is a constant, too, similarly to the case of two metallic 
spheres. The peculiarity of our case is that the presence 
of monovalent salt which screens the surface electric po- 
tential. Deep inside r s -cylinder, screening by monovalent 
salt is incomplete. Indeed, according to the Debye-Huckel 



— —J - 

^ 





' / r/ 




7\ 


2W S 




j 




< 








" — ^ 
— >- 


2r s r<— 



FIG. 7: The cross section of the two aggregated spherical 
macroions. Z-ions adsorbed on them are not shown. The 
width of the gray region is the effective range of the Coulomb 
interaction in the presence of monovalent salt. The neutral rs- 
cylinder is shown by the rectangle of dashed lines. It is larger 
than the sticky region shown by the full line rectangular (see 
Fig.©. 



theory, the density of screening charge of monovalent 
salt near the macroion surface has the same magnitude 
p = —ip/Anr^ as outside of r s -cylinder. But with decreas- 
ing d, the volume occupied by this charge shrinks very 
fast. At d <C r s , charge density of monovalent salt pro- 
jected to the macroion surface, a' — —(fid/Airr^ is much 
smaller than <7o . As we saw from the problem of two con- 
tacting metallic spheres, the total surface charge density 
between two touching spherical surfaces with constant 
potential decays very fast towards the contact point. Ac- 
cordingly, inside the r s -cylinder at d <C r s , the surface 
charge density a + cr'/2 is very small (it decays with de- 
creasing d as <7o exp(n(y l / R/2r s — y/R/d))). Therefore, 
for the surface charge density of the macroion a{d), we 
arrive at aid) = —a'(d)/2 = 4>d/8irr^. Since at d <C r s , 
u(d) is much smaller than cro, we call this region "neu- 
tral region" . At r <C r s , the sticky region discussed above 
(the full line rectangular in Figs. HH and 01 is much smaller 
than the neutral region and therefore almost neutral. 

Let TV be the number of Z-ions in a free macroion- Z- 
ions complex. We define Q* = -Q + NZe = Ze(N — 
Ni) as the net charge of a free complex, and (1 — fi)Q* 
as the net charge of the aggregated complex. The total 
number of Z-ions in an aggregated complex is therefore 
N — (3(N — Ni). 

At r s ^> R, the whole complex is neutral in the aggre- 
gate, so (3 = 1. At r s <§C R, away from neutral regions, 
due to equilibration of Z-ions, aggregated macroion-Z- 
ions complexes have the same surface density of Z-ions 
as free complexes. Thus, (3 is a measure of the fraction 
of total surface area of a complex occupied by neutral 
regions. j3 decreases with decreasing r s . 

Now we can write down the free energy of the system. 
Similarly to Ref. GH we start from a point (p, s) in the 
domain of partial aggregation and assume that a fraction 
x of aggregated macroion- Z-ions complexes. The free 
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energy per unit volume is 



where 



F (N,x;s,p) = (l-x)p 



O* 2 

^+NE(N) 



+ xp 



(1 - 0)^ + NE(N) - Pfi c (N -N,) + e 



+ k B T[s - (1 - x)pN - xp(N - (3(N - Ni))] 

. [a - (1 - x) P N - xp{N - (3{N ~ Nj))]vo 
x in . 



(i) 



The first line of Eq. JTJ is the free energy density of free 
macroion- Z-ions complexes. It consists of Coulomb en- 
ergy of complexes and correlation energy of Z-ions. Here 
C is the capacitance of a free complex, E(N) < is 
the correlation energy per Z-ion in the complex. The 
second line of Eq. is the free energy density of ag- 
gregated complexes. Here e < is the additional cor- 
relation energy per complex due to aggregation^ and 
fx c (N) = d(NE(N))/8N < is the part of the chemical 
potential of Z-ions related to their correlations on the 
surface of a macroion. It has been assumed in Eq. (JTJ 
that N — Ni <C Ni (this will be confirmed later) and 
therefore 



Hc(N) = n c {Ni) = n c 



(2) 



In the second line of Eq. the first term is the Coulomb 
energy of an aggregated complex diminished by a loss of 
fraction (3 of capacitor area, the second and third term 
give correlation energy of Z-ions diminished by loss of 
j3(N — Ni) of them and the fourth term is responsible 
for attraction energy of aggregated complexes in sticky 
regions. The third and fourth line of Eq. Q give the free 
energy density of free Z-ions related to their entropy, (vq 
is the normalizing volume). The entropy of macroions 
is ignored because of their much smaller concentration. 
Minimizing this free energy with respect to N and x, we 
get 



ZeQ* 
C 



k B T\a[{s - (l-x)pN 
xp{N - (3(N - Ni)))v ], 



(3) 



e-/3^ = 0(N-N i ){/ J i c -kBT\n[(s-(l-x)pN 



xp(N - 0(N - Ni)))v }}, 



(4) 



Excluding the entropy term from in Eqs. J3J| and JIJ, 

we get 



N Ctd = Nil It 



2C(r s ) 



0(r s ) Q 2 



(5) 



Here and below, the upper (lower) sign in the formula 
always corresponds to the first (second) subscript of the 
symbol. 

Plugging Eq. JSJ back to Eq. J3J), we get 

s c ,d(p) = s c , d +(l-x)pN Cid +xp[N C}d -0(r s )(N c , d -N i )], 

(6) 



Sc,d 



t'o 



■ exp 




(7) 



The two solutions N c<d and s c ^ d mean that we have two 
partial aggregation domains. Taking x = and x = 1 in 
Eq. |JBJ, we get two outer boundaries 



s c ,d(p; x = 0) = s cA + pN c>d , 
and two inner boundaries 



(8) 



s c , d (p; x = 1) = s c . d + p[N c<d - 0{r s ){N c . d - Ni)] (9) 

of the partial aggregation domains. The corresponding 
phase diagram is shown in Fig. 2t>. The outer bound- 
aries are also the boundaries of the whole aggregation do- 
main (including partial and total aggregation domains), 
at which all free macroion- Z-ions complexes are still sta- 
ble (a; = 0). The inner boundaries are also the bound- 
aries of the total aggregation domain at which there is 
no free macroion- Z-ions complex (x = 1). When s in- 
creases from zero, first the aggregate forms, then it dis- 
solves. Thus we arrive at a phase diagram of reentrant 
condensation. According to Eq. @, the slopes, ds/dp, 
of the outer boundaries are N Ctd , which are the num- 
bers of Z-ions in a free macroion- Z-ions complex. Since 
N Cjd are solutions to N in the whole partial aggregation 
domain, the number of Z-ions complexed with one free 
macroion is fixed in the partial aggregation domain. Ac- 
cording to Eq. ©, the slopes of the inner boundaries 
are N Ctd — 0(N c>d — Ni), which are the numbers of Z- 
ions in an aggregated complex. On the other hand, s Ctd 
are the intercepts of these boundaries with p axis. The 
dashed line on Fig.^J) corresponds to solutions with neu- 
tral macroion- if-ions complexes and can be calculated 
from Eq. J3J| . Taking N = Ni in this equation, we get 



i ( lMd\ 

= — cxp - — 



so{p) = — exp ( I +pNi. 

Vq \ k B i I 



(10) 



This line has the same slope as the isoelectric line Si = 
pNi (dotted line in Fig. but a small finite intercept. 
This is because there is always a small fraction of free 
Z-ions in the solution. 

When 0=1, the whole surface of a macroion in the 
aggregate is neutral, all formulae above reproduce the 
corresponding results in the case r s i^ii. For example, 
the boundaries of the total aggregation domain given by 
Eq. are now parallel lines 



s c ,d(p; x = l) = s c , d + pNi. 



(11) 



The corresponding phase diagram for r s 3> R is shown 
in Fig. 

There are four parameters /i c , e, C and (3 in the theory. 
Their values depend on the screening radius r s , the shape 
and flexibility of the macroions and Z-ions. In the follow- 
ing sections, we calculate them for specific systems. But 
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even now we can qualitatively summarize the evolution of 
the phase diagram with decreasing r s . For this purpose, 
we define the size of an aggregation domain as the abso- 
lute value of the difference of the slopes of its two bound- 
aries. As r s decreases, (3{r s ) decreases, C(r s ) increases, 
while e is fixed (this will be shown below). Then accord- 
ing to Eq. (J5J, Nd increases and iV c decreases. Thus, 
both the size of the whole aggregation domain (including 
partial and total aggregation domains), Nd — N c , and the 
size of the total aggregation domain, (1 — (3){N d — N c ), 
grow with decreasing r s . In contrary to them, the rel- 
ative size of the two partial aggregation domains in the 
whole aggregation domain, (3{Nd — N c )/(Nd — N c ) = (3, 
decreases with decreasing r s . 

As we mentioned before, at r s < r the growth of the 
aggregation domain and reduction of the relative size of 
the partial aggregation domain continue with decreasing 
r s . In this regime, however, the range of Coulomb in- 
teraction and correlation becomes the same and one can 
not separate these interactions from each other. A micro- 
scopic theory is necessary, which operates with screened 
interaction of Z-ions with macroion and with each other 
and allows for small dielectric constant of macroion. This 
is a subject of future work. 



III. RIGID SPHERICAL MACROIONS 

In this section, we consider the system with spherical 
macroions and much smaller spherical Z-ions (Figs. 
We estimate parameters fi c , e, C and (3 microscopically 
and find out how borders of the aggregation domain 
change with decreasing screening radius r s . 

Let us first consider the case r s R. In this case, 
each macroion carry iV, Z-ions in the aggregate and 
(3=1. In the first order approximation, Z-ions repel 
each other to form 2D Wigner crystal on the surface of a 
macroion. Thus the correlation chemical potential /i c is 
approximately^ 



fj, c 



1.6ZV 
Dr 



(12) 



where r is the radius of the Wigner-Seitz cell when a 
macroion is neutralized by Z-ions (see Figs. |5J and D is 
the dielectric constant of water. To calculate e, we first 
notice that the number of Z-ions in a sticky region is 



Ni 



nW 2 
4vri? 2 



(13) 



where relation Ni = 4ttR 2 /irr 2 has been used. If each 
macroion has 6 nearest neighbors in the aggregate (6=12 
for dense packing), 



'Nidfjic = —0.467a 



ZeQ 
DR ' 



(14) 



Here we assumed that each Z-ion in the sticky region 
gains the correlation energy a(i c . The maximum possible 



a can be estimated as follows. Since the surface density 
of Z-ions is doubled in the contact region, the radius of 
the Wigner-Seitz cell becomes and the correlation 

energy per Z-ion gets a factor y/2. This gives a ~ 0.3. 
The real a is much smaller than this upper limit due to 
the effect of small dielectric constant of macroions. As 
is well known in electrostatics, positive Z-ions in water 
(D ~ 80) creates positive images in adjacent macroions 
(D <C 80). Repulsion from images push two contacting 
macroions away from each other and diminishes the gain 
of correlational energy. Finally, for a spherical macroion 
with Z-ions whose size is negligible, the capacitance is 



C = DR. 
Using Eqs. © and 0, we get 

N c , d = Ni (l T V0-26 7 a^ 

1.6ZV 
k B TDr 



Sc,d 



— exp 



1± 



67a 
"08" 



(15) 

(16) 
(17) 



Let us now switch to the case r -C r s <C R. Since the 
sticky region is much smaller than the neutral region (see 
Fig. EJ , the short-range correlations are not changed by 
screening, fi c and e are still given by Eqs. iff^l and fH)) . 
But (3 and C are affected by screening. For simplicity, 
we estimate (3 by assuming that the whole r s -cylinder is 
neutral. Then j3 reduces to the area ratio of the neutral 
region to the whole surface of the complex, 



(3 = b 



br„ 



AttR 2 2R' 



And 



This gives 



C = 



DR 2 



N c , d = Nil It ^0.4 7 a- ) , 



— exp 




k B TDr 



(18) 
(19) 

(20) 
(21) 



It is easy to check that in both cases N C}d are close to 
Ni. This justifies the approximation J5J used in the last 
section. From Eq. H20(l . we see that the aggregation do- 
main broadens with decreasing r s and the angle it oc- 
cupies becomes the order of unity at r s ~ r. According 
to Eq. (fT%|) . the relative size of the partial aggregation 
domain decreases with r s proportionally to r s /R. 

Two comments on validity of our approach are in order 
here. First, we assumed above that at s > s d (p), when 
overcharged macroions redissolve, they still keep almost 
all Z-ions. As is well knownSi, the condition of maximum 
charge inversion is given by 



ZeQ* 
C 



(22) 
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Denoting N max the maximum number of Z-ions on a 
macroion, and using Eq. (|12[l. in the case of r s <C R 
we have 



AT = Ni 1 + 0.4 



(23) 



Comparing with Eq. I|2U|) , we see that N max and Nd are 
of the same order of magnitude. As we mentioned above 
we believe that a < 0.3 so that Nd < N max and the 
picture of stable complexes on both sides of aggregation 
domain is justified. 

Second, up to now we discussed screening by mono- 
valent salt in Debye-Hiickel approximation. When a Z- 
ion has large charge and small size a, they are screened 
nonlinearly. Monovalent counterions condense on Z-ions 
(similarly to the Manning screening of a charged cylin- 
der) and renormalize the charge of a free Z-ion to a 
smaller value Z— . The charge Z' is screened in Debye- 
Hiickel way. When adsorbed to the macroion surface, Z- 
ion is subjected to the additional Coulomb interactions 
with both the negative macroion surface and neighbor- 
ing Z-ions. Both interaction energies are of the order of 
Z'e 2 /r and at a <C r can be ignored in comparison with 
with the chemical potential of the monovalent counterion 
Z'e 2 /a on Z-ion. Therefore, the effective charge Z' is not 
changed by Z-ion condensation and should be used as the 
effective charge of a Z-ion in all our results 2 ^. 

In recent Monte Carlo simulations of the system of 
spherical macroions with multivalent Z-ions 24 , it is seen 
that around the isoelectric line Sj = pNi, all macroions 
aggregate, while far away from this line on both sides, free 
macroion- Z-ions complexes are stable. This qualitatively 
agrees with our theory. 



IV. CYLINDRICAL MACROIONS 

In this section, we consider the system with cylindrical 
macroions and much smaller spherical Z-ions (Fig. I3J. 
We show that the qualitative feature of the phase dia- 
gram of this system is the same as spherical macroions. 
The major difference is that for cylinders, Sd is much 
larger. 

We assume that charge of the macroion and Z-ion is 
such that r <C R. Accordingly two dimensional Wigner- 
Crystal-like structure of Z-ions can be formed on the 
surface of macroions 19 . We also assume that r s < L. 
Consequently below we focus on two cases corresponding 
to R <C r s <C L and r « r s « Also note that in order 
to gain more contact area and, therefore, more correla- 
tion energy, all cylinders in the aggregate are parallel to 
each other. 

We first consider the case R <C r s <C L. It is easy 
to see that the aggregate should be neutral, i.e., (3=1, 
and the chemical potential /i c is still given by Eq. (|12fl . 
However, e is not the same since the area of the sticky 
region is much larger. It is not a disk, but a stripe with 



area 2WL. The number of Z-ions in the sticky region is 
2WL 



_ . . _ 1 ryr 

ZnRL 7T V R 



(24) 



For a macroion with b nearest neighbors (6 = 6 for dense 
packing) , the additional correlation energy per macroion 
is 



ba fyr 
e = — NiiJ, c * —. 

TT V R 



(25) 



Here a takes the same value as the one for spherical 
macroions. And the capacitance is 



C 



DL 



2ln(r s /R)- 
Using Eqs. JSJ and (J7J), we have 



Nr d = Ni 



(26) 



3rd 



— exp 




(27) 



In^f* 
0.2tt R V r 



(28) 



Now we consider the case r <§; r s <§; R. Instead of a 
neutral aggregate, we only require a stripe-like neutral 
region in the present case. It is easy to see that both /z c 
and e are not changed by screening. The fraction of the 
neutral area is 



p = b 



2W,L b 2r« 



2-kRL tt V R 



while the capacitance 



C 



DRL 

2r, 



We have 
N c .d = Ni 



1=F \ 0.8a 



7N 2 



1 

— exp 



2 

1.6Z 2 e 2 
k B TDr 



1 ± 



0.2 V2/ ' r 



(29) 



(30) 



(31) 



(32) 



Again in both cases N Ci d are close to iVj and the ap- 
proximation J2J is valid. According to Eq. (|29(l , the rela- 
tive size of the partial aggregation domain decreases with 
r s proportionally to ^r s /R. We see that the aggregation 
domain are wider in the case of cylindrical macroions 
than that for spherical macroions (see Sec. III). This 
happens because of the larger contact area in the former 
case. 
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According to Eq. (|32|l . formally Sd can be even larger 
than 1/vq. Consequently, one may conclude that the ag- 
gregate never dissolve at large s. Actually, this is not 
true. When s becomes so large that the distance be- 
tween Z-ions in the bulk is equal to r s , the repulsion 
between free Z-ions can not be ignored and our theory 
is invalid. At such a large concentration, Z-ions are cor- 
related not only on the surface of macroions but also in 
the bulk and aggregation does not happen. What we can 
conclude is that for systems discussed in this section, the 
aggregate dissolves at very large s. Mathematically, s d 
is so large that ln(l/sd«o) < ln(sd/sc)- Correspondingly, 
the aggregation domain in Fig. 0] at s > is very wide. 

The results above can be qualitatively compared with 
experimental results for solutions of DNA with spermine 
in which long DNA double helices may be consid- 
ered as rigid negative cylinders and short positive sper- 
mine molecules as small Z-ions (Z — 4). In experiments, 
reentrant condensation is observed. Also the aggregation 
domain at small p is very wide (log(sd/s c ) ~ 4). Both 
facts agree with our theory. 



V. ARTIFICIAL CHROMATIN 

In this section, we discuss the artificial chromatin in 
which macroions are long semi-flexible polymers with lin- 
ear charge density — rj and Z-ions are hard spheres with 
radius R (Fig.^l. We assume Rrj <C Ze, i.e., many turns 
of polymer segments are needed to neutralize one sphere. 
The aggregate of such a system looks like Fig. 3 in Ref.lTll 
The phase diagram of this system has already been dis- 
cussed in detail in Ref. ^] under the assumption that 
the aggregate is always neutral which is valid only for 
r s ^> R. When r s <C R, there is no justification for 
the neutrality of the aggregate and the theory in Sec. ITT1 
should be used. Here we calculate the parameters mi- 
croscopically and obtain the phase diagram for the later 
case. 

First, let us remind the results following directly from 
Ref. 1 1 1 L In this case, each macroion- Z-ions complex car- 
ries Ni Z-ions in the aggregate. Correspondingly, (3=1. 
The correlation chemical potential [i c is essentially the 
self-energy of a bare free sphere in the solution which is 
almost totally eliminated in the complex. Therefore 



2^2 



Z 2 e 
2DR ' 



(33) 



where D is the dielectric constant in water solution. To 
calculate e, we first notice that the correlation energy 
per unit length of the polymer is just the interaction en- 
ergy of the polymer turn with its stripe of background 
(sphere) positive charge. Thus, it is — rj 2 ln(R/A)/D. 
Correspondingly, when the density doubles, A is halved, 
the gain in correlation energy per unit length is —ar/ 2 /D, 
where a = In 2 as a rough estimate. Similar to the case 
of spherical macroion, the radius of the sticky region 
W ~ y/RA/2. And the polymer length in this region 



is rrW 2 /A = ttR/2. Since each aggregated complex car- 
ries Ni Z-ions and has b nearest neighbors, we have 



irba Rr/ 2 
~2 ~D 



(34) 



where b = 6 if the maximum packing number is achieved. 
The capacitance 



C 



DL 



DRN,; 



2\n{r s /R) ln(r s /i?)' 



(35) 



where we have used L ~ 2RNi, because near neutrality, 
each polymer absorbs about Ni spheres and all spheres in 
the macroion- Z-ions complex are densely packed. Using 
Eqs. (J5J and J7J), we have 



Nr. 



N, 
1 



it 



■ exp 



/ irba Rrj 
\n{r s /R)~Ze 

Z 2 e 2 



(36) 



. I±2j7r6aln^^ 
2k B TDR V V RZe 



(37) 



This gives a typical phase diagram of reentrant conden- 
sation (Fig.^Ji). 




2W. 



FIG. 8: A 2D view of the surface of a sphere touching the 
other sphere at the center of the disk-like region a. To make 
a neutral, a semi-flexible polymer have to change the density 
of its turns to make the whole equatorial stripe-like region b 
neutral. The rest of the surface colored in gray is charged. 

Now let us consider the case when A <C r s -C R- A 
being the average distance between two polymer turns 
on the surface of spheres when polymer segments exactly 
neutralize the sphere (see Fig. QJ. Since rj/A ~ Ze/R 2 , 
the requirement A <C R is equivalent to Rrj <C Ze. As 
we discussed in Sec. ^ with the screening of monova- 
lent salt, the whole aggregate need not be neutral, but 
disk-like regions on spheres where two spheres touch each 
other must be neutral (see Fig. UJ. The radius of this 
region W s = y/2r s R. In the present case, the neutral 
region should be made by the rearrangement of the poly- 
mer segments on the surface of spheres. If the polymer 
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is extremely flexible, it can bend easily in the disk-like 
region. By this bending, the surface density of polymer 
segments is changed "locally" so that the region can be 
made neutral and the rest of the surface is not changed. 
However, if the polymer is only semi-flexible (the persis- 
tent length l p is such that \/r s R <C l p <C R) , it is too rigid 
to neutralize the disk-like region. Then this region can be 
made neutral only by changing distances between sequen- 
tial turns of the polymer on the surface of the sphere. For 
example, for two touching spheres (Fig.0, the equatorial 
stripe-like region shown in Fig. [S] can be made neutral. 
We focus on the case of the semi-flexible polymer which 
corresponds to DNA. 

When A <C r s -C R, the correlation chemical potential 
fi c is still given by the self-energy of a bare free sphere, 
which is now 



/'c 



r s Z 2 e 2 
2DR 2 ' 



(38) 



and the correlation gain due to aggregation, e, is still 
given by Eq. (|34|) . This is because the contact region 
where the additional correlational energy is gained is 
much smaller than the neutral region as discussed in 
Sec. [n] Also since the number of U-ions in each com- 
plex in the aggregate is close to N, in the first order 
approximation, we can use Ni as the number of Z-ioms 
in each complexes in Eq. I|34fl ■ Knowing that the area of 
the strip- like region is 4irW s R (see Fig.JSJ, the area ratio 
(3 is given by 



4ttW s R 
AttR 2 



And the capacitance is now 

c= drl = d En ^ 

2r s r s 
Using Eqs. JSJ and Q, we have 



(39) 



(40) 



than the correlation energy of the sphere in the complex 
(it interacts with all turns of the polymer) . 

An artificial chromatin system in the presence of mono- 
valent salt is studied in experiments^, where the phase 
diagram of reentrant condensation is found and the do- 
main of partial aggregation is not observed. This is qual- 
itatively consistent with our theory that the relative size 
of the partial aggregation domain is small due to screen- 
ing effect of monovalent salt. 



sZe 



= Ni 



it 



i 

— exp 



irbaRrj (R 

r s Z 2 e 2 
2k B TDR 2 



1±2 4 



' irba Rrj f R 



(41) 



(42) 



The phase diagram is shown in Fig. 0b . 

It is easy to check that in both two cases, slopes N Cf d 
are close to N and the approximation is valid. Com- 
paring s c .d and N c ,d with results for spherical and cylin- 
drical macroions, we see that aggregation domain of the 
artificial chromatin is the smallest of all three cases if pa- 
rameter Rrj/ Ze is small (the polymer makes many turns 
around the sphere). The main reason is that in artifi- 
cial cromatin, e and fi c originate from two different kinds 
of correlations. At Rr/ <C Ze, the correlation energy of 
two polymer turns of the touching sphere is much smaller 




FIG. 9: The chromatin regulation and compaction. sZe is the 
concentration of charge of histones, pQ is the concentration 
of DNA. All lines and charges have the same meaning as in 
Fig. 0] Normal chromatin is at the point 1. Acetylation of 
histones, moves it to the point 2 (downward arrow). On the 
other hand, addition of protamines during spermatogenesis 
compacts chromatin into the more condensed state 3 (upward 
arrow) . 



The phase diagram Fig. 0Jd may be useful for qualita- 
tive understanding of chromatin regulation. It is known 1 
that chromatin can switch from the compact state of 30 
nm fiber or higher order structures to the loose state 
of 10 nm fiber in which gene transcription takes place. 
This transition is caused by acetylation of core histones, 
which reduces the number of positive charges of histone 
octamers. Experiments in vitroi^ show that the transi- 
tion happens when the charge of the octamer Ze (about 
+160e in normal conditions) is reduced by +12e, i.e., 
only by 8%. Such a sensitive response can be understood 
in our theory as shown in the phase diagram Fig. [§] in 
the plane of total charges of octamers and DNA. Nor- 
mal chromatin seems to be partially aggregated and the 
system is at state 1. Under acetylation, the charge of 
octamers is reduced so that the system moves down to 
state 2 below the partial aggregation domain. The small 
partial aggregation domain is crucial for the sensitive re- 
sponse. Fig. [§] also helps to understand spermatogene- 
sis in which normal chromatin is compacted further by 
adding strongly positively charged proteins, protamines 1 . 
In the solution, the concentration of the positive charge 
increases and the system moves up from state 1 to state 
3 in the total aggregation domain. 
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To conclude this section, let us make a comment about 
the role of Manning condensation of monovalent salt on 
DNA. Above we literally dealt with the case when the 
linear charge density of the polymer rj is smaller than 
critical density r/ c — ksT/e of Manning condensation. 
As is well known, for a free polymer with rj > rj c (for 
example, for DNA double helix, which has 77 = 4.2r/ c ), 
X] is renormalized to r\ c . For DNA wrapping a positive 
sphere in artificial chromatin, this renormalization is dif- 
ferent because of the surface charge of positive spheres. 
Such renormalization and its consequence on the phase 
diagram were studied in Ref. [nl 

VI. CONCLUSION 

In this paper we developed a phenomenological the- 
ory of complexation and aggregation in asymmetric wa- 
ter solutions of strongly charged negative and smaller 
positive colloids (macroions and Z-ions) in the presence 
of a large concentration of monovalent salt. We showed 
that in contrary with earlier theoryii mostly devised for 
salt free solutions, aggregate of complexes can carry net 
charge of macroions and Z-ions which is almost as large 



as the charge of free complexes. Only small touching re- 
gions of each aggregated complex are depleted of the net 
charge. As a result screening by monovalent salt leads to 
broader aggregation domain in the phase diagram of the 
solution and to narrower, but still finite domain of par- 
tial aggregation where aggregate is in equilibrium with 
free complexes. Our phenomenological theory expressed 
properties of the phase diagram through several param- 
eters which depend on charges of macroions and Z-ions, 
their sizes, shapes, flexibility and screening of monovalent 
salt. For three different pairs of macroions and Z-ions we 
evaluated these parameters and discussed differences and 
trends in their phase diagrams. We found out that one 
of these systems, which we call artificial chromatin (see 
Fig- Q) j indeed has phase diagram qualitatively similar 
to the one of natural chromatin. 
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